arXiv:1507.00359v2 [hep-ex] 29 Dec 2015 


EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2013-037 

2015/12/31 


CMS-HIG-14-025 


Search for exotic decays of a Higgs boson into undetectable 
particles and one or more photons 


The CMS Collaboration^ 


Abstract 

A search is presented for exotic decays of a Higgs boson into undetectable particles 
and one or two isolated photons in pp collisions at a center-of-mass energy of 8 TeV. 
The data correspond to an integrated luminosity of up to 19.4 fb -1 collected with the 
CMS detector at the LHC. Higgs bosons produced in gluon-gluon fusion and in as¬ 
sociation with a Z boson are investigated, using models in which the Higgs boson 
decays into a gravitino and a neutralino or a pair of neutralinos, followed by the de¬ 
cay of the neutralino to a gravitino and a photon. The selected events are consistent 
with the background-only hypothesis, and limits are placed on the product of cross 
sections and branching fractions. Assuming a standard model Higgs boson produc¬ 
tion cross-section, a 95% confidence level upper limit is set on the branching fraction 
of a 125 GeV Higgs boson decaying into undetectable particles and one or two isolated 
photons as a function of the neutralino mass. For this class of models and neutralino 
masses from 1 to 120 GeV an upper limit in the range of 7 to 13% is obtained. Further 
results are given as a function of the neutralino lifetime, and also for a range of Higgs 
boson masses. 
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1 Introduction 

The detailed studies of the properties of the observed Higgs boson |Tj-j3l are key components 
of the LHC physics program. In the standard model (SM) and for a given mass of the Higgs 
boson, all properties of the Higgs boson are predicted. Physics beyond the SM (BSM) might 
lead to deviations from these predictions. Thus far, measurements of the Higgs bosons cou¬ 
plings to fermions and bosons and of the tensor structure of the Higgs boson interaction with 
electroweak gauge bosons show no significant deviations 00 with respect to SM expectations. 

Measurements of Higgs boson couplings performed for visible decay modes provide con¬ 
straints on partial decay widths of the Higgs boson to BSM particles. Assuming that the cou¬ 
plings of the Higgs boson to W and Z bosons are smaller than the SM values, this indirect 
method provides an upper limit on the branching fraction of the 125 GeV Higgs boson to BSM 
particles of 57% at a 95% confidence level (CL) 00. An explicit search for BSM Higgs boson 
decays presents an alternative opportunity for the discovery of BSM physics. The observation 
of a sizable decay branching fraction of the Higgs boson to undetected (e.g. invisible or largely 
invisible) final states would be a clear sign of BSM physics and could provide a window on 
dark matter |7-10|. 

Several BSM models predict Higgs boson decays to undetectable particles and photons. In 
certain low-scale supersymmetry (SUSY) models, the Higgs bosons are allowed to decay into a 
gravitino (G) and a neutralino (x?) or a pair of neutralinos jlTl 12 ]. The neutralino then decays 
into a photon and a gravitino, the lightest supersymmetric particle and dark matter candidate. 
Figure [l] shows Feynman diagrams for such decay chains of the Higgs boson (H) produced by 
gluon-gluon fusion (ggH) or in association with a Z boson decaying to charged leptons (ZH). 



Figure 1: Feynman diagrams for the H —> undetectable + 7 final state produced via ggH (left) 
and ZH (right). 

As the gravitino in these models has a negligible mass ITTI 121. the remaining parameter is the 
neutralino mass. If its mass is in the range ?»h/2 < 0 < ni\\, with mu = 125 GeV the mass 

of the observed Higgs boson, the branching fraction £>(H —>■ x^G —» 7 GG) can be large. For 
m -i) < m H /2, the decay H —> x^Xi H► 77 GG is expected to dominate. The same discussion can 
be applied to heavy neutral Higgs bosons with masses larger than 125 GeV. The lifetime of the 
neutralino can be finite in some classes of BSM scenarios, leading to the production of one or 
more photons displaced from the primary interaction. 

In the SM, the signature equivalent to the signal arises when the Higgs boson decays as H —> 
Z 7 —> i/i /7 with a branching fraction of 3 x 10 4 . The decay H —> Z 7 has been studied in 
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2 The CMS experiment 


Z —> e 1 e and Z —> pi + }i~ final states. Upper limits on the product of the cross section and 
branching fraction of about a factor of ten larger than the SM expectation have been set at the 
95% CL IlSillJ- With the available dataset the search presented is not sensitive to this decay, 
but it is sensitive to enhancements in the Higgs boson decay rates to undetectable particles and 
photons arising from BSM physics. 

Various background processes lead to the signal signatures and are estimated from simulation 
or from control samples in data. The dominant background processes are from y+jets events 
and diboson events in the ggH and ZH search, respectively. Details of the background estima¬ 
tion techniques are discussed in Section[5] The strength of the ZH channel analysis is an almost 
background-free selection leading to a larger sensitivity in the model-dependent interpretation. 
While both the ggH and the ZH channels provide sensitivity to BSM Higgs boson signatures, 
the ggH channel allows a model-independent interpretation of the results. 

This analysis presents a first search for decays of a scalar boson to undetectable particles and 
one or two isolated photons. The scalar boson is produced in ggH or in ZH. The data used 
correspond to an integrated luminosity of up to 19.4 ± 0.5 fb 1 at a center-of-mass energy of 
y/s = 8 TeV in 2012 collected with the CMS detector at the CERN LHC. 

The results of the search are presented in terms of the low-scale SUSY breaking model for 
» 4 -[ = 125 GeV and m-u between 1 GeV and 120 GeV, and for mu between 125 GeV and 400 GeV 

Ai 

for the example case where m^o = mu — 30 GeV. The effect of a finite lifetime (t~o) is studied 

for the example case where 7 «h = 125 GeV and m~a = 95 GeV. 


2 The CMS experiment 

The CMS detector, definitions of angular and spatial coordinates, and its performance can be 
found in Ref. Ifl5ll . The central feature of the CMS apparatus is a superconducting solenoid, of 
6 m internal diameter, providing a magnetic field of 3.8 T. The field volume contains a silicon 
pixel and strip tracker, a crystal electromagnetic calorimeter (ECAL), and a brass and scintilla¬ 
tor hadron calorimeter. Muons are measured in gas-ionization detectors embedded in the steel 
flux-return yoke of the magnet. The first level of the CMS trigger system, composed of special¬ 
ized hardware processors, is designed to select the most interesting events within 3 pis, using 
information from the calorimeters and muon detectors. A high-level trigger processor farm is 
used to reduce the rate to a few hundred events per second before data storage. 

A particle-flow algorithm [16j, 17] is used to reconstruct all observable particles in the event. 
The algorithm combines all subdetector information to reconstruct individual particles and 
identify them as charged hadrons, neutral hadrons, photons, and leptons. The missing trans¬ 
verse energy vector E™ ss is defined as the negative vector sum of the transverse momenta of 
all reconstructed particles (charged or neutral) in the event, with Ej' ss = \Ej ' ss \. Jets are re¬ 
constructed using the anti-kj clustering algorithm HT8l with a distance parameter of R = 0.5, as 
implemented in the FASTjET package l!T9l 12011 . A multivariate selection is applied to separate 
jets from the primary interaction and those reconstructed due to energy deposits associated 
with pileup interactions |2TI . The discrimination is based on the differences in the jet shapes, 
on the relative multiplicity of charged and neutral components, and on the different fraction 
of transverse momentum which is carried by the hardest components. Photon identification 
requirements and other procedures used in selecting events can be found in Section [4] 
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3 Data and simulation events 

In the search for Higgs bosons produced in ggH, the trigger system requires the presence of 
one high transverse energy (Ej) photon candidate and significant E™ 1SS . The presence of a pho¬ 
ton candidate with Ej > 30 GeV is required within the ECAL barrel region (\if f \ < 1.44). At the 
trigger level E™ ISS is calculated from calorimeter information, and is not corrected for muons. 
A selection requirement of E™ 1SS > 25 GeV is applied. The efficiency of the trigger is monitored 
and measured with two control triggers for the photon and the E™ 1SS trigger requirement. The 
data recorded with this trigger correspond to an integrated luminosity of 7.4 fb 1 and were 
part of the CMS "data parking" program implemented for the last part of the data taking at 
y/s = 8 TeV in 2012. In that program, CMS recorded additional data with relaxed trigger re¬ 
quirements, planning for a delayed offline reconstruction in 2013 after the completion of the 
LHC Run 1. 

For the search for Higgs bosons produced in ZH, collision events were collected using single¬ 
electron and single-muon triggers which require the presence of an isolated lepton with pj in 
excess of 27 GeV and 24 GeV, respectively. Also a dilepton trigger was used, requiring two 
leptons with pj thresholds of 17 GeV and 8 GeV. The luminosity integrated with these triggers 
at \Js = 8 TeV is 19.4 ftV 1 . 

Several Monte Carlo (MC) event generators are used to simulate signal and background pro¬ 
cesses. The simulated samples are used to optimize the event selection, evaluate selection 
efficiencies and systematic uncertainties, and compute expected event yields. In all cases the 
MC samples are reweighted to match the trigger efficiency measured in data. 

The V 7 , WZ, ZZ, VW (where V represents W or Z bosons), Drell-Yan (DY) production of qq —> 
Z/ 7 *, and qq —» W 1 W processes are generated with the MadGraph 5.1 event generator 1221 
at leading-order (LO), the gg —> W 1 W ~ process is generated with the LO event generator 
GG2WW 3.1 [ f23l , and the tt and tW processes are generated with POWHEG 1.0 at next-to-leading- 
order (NLO). The signal samples are also produced with MadGraph. The cross sections at 
NLO or higher orders if available are used for a given process to renormalise the MC event 
generators. All processes are interfaced to the PYTHIA 6.4 generator Il24l for parton shower and 
hadronization. 

The CTEQ 6 L set of parton distribution functions (PDF) Il25l is used for LO generators, while 
the CT10 H26H PDF set is used for NLO generators. For all processes, the detector response is 
simulated with a detailed description of the CMS detector, based on the Geant4 package H2 7\ . 
Additional pp interactions overlapping the event of interest in data, denoted as pileup events, 
are accounted for by simulating pp interactions with the PYTHIA generator and adding them 
to each MC sample. The MC samples are tuned to reproduce the distribution in the number of 
pileup events in data. The average number of pileup events is about 26 for the collected data 
used in the ggH channel, and is about 21 for the collected data used in the ZH channel. 


4 Event selection 

Two strategies are followed to isolate the Higgs boson events produced by ggH and by ZH from 
the background processes. The signal cross sections are several orders of magnitude smaller 
than the major reducible background processes, whose contributions are greatly reduced using 
the event selections described in the following sub-sections. 
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4 Event selection 


4.1 Event selection in the ggH channel 

In the ggH channel, each selected event is required to have at least one photon candidate with 
Ej > 45GeV and \i]' r \ < 1.44 using a cut-based selection H28ll29l . To reduce the SM back¬ 
grounds arising from the leptonic decays of W and Z bosons, a lepton veto is applied. Events 
are rejected if they have one or more electrons fulfilling a loose identification requirement l30l 
and Pj > lOGeV and |?/ e | < 2.5, excluding the transition region of 1.44 < \p e \ < 1.57 since the 
reconstruction of an electron object in this region is not optimal. Similarly, events containing 
muon candidates with pj > lOGeV, \if‘\ < 2.1, and well separated from the photon candidate 
requiring AR(y, p) = V(A //) 2 + (A <p) 2 > 0.3 (where the cp is azimuthal angle in radians), are 
rejected. In addition to the selection requirements described above, the E™ ISS is required to be 
greater than 40 GeV. This level of selection is referred to as the preselection. Additional selec¬ 
tion criteria are applied to search for new physics in either a quasi model-independent way or 
optimized for a SUSY benchmark model. In this channel jets can arise from initial-state radia¬ 
tion. For both search strategies jets are required to have Pj > 30 GeV and rji \ < 2.4. These jets 
must not overlap with the photon candidate below AR( 7 ,jet) < 0.5. 

In the model-independent analysis, events with two or more jets are rejected. For events with 
one jet the azimuthal angle between the photon and the jet (A<^>('y, jet)) is required to be smaller 
than 2.5. This selection requirement rejects the dominant 7 +jet background, where the photon 
and the jet tend to be back to back in the transverse plane. 

In the model-dependent analysis developed for SUSY scenarios, no requirement is applied 
on jet multiplicity. In order to minimize the contribution from processes such as 7 +jets and 
multijet events, two methods are used for identifying events with mismeasured Ey ISS . The 
Ey ISS significance method II3T1 takes account of reconstructed objects for each event and their 
known resolutions to compute an event-by-event estimate of the likelihood that the observed 
Ey ISS is consistent with zero. In addition, a minimization method |[29l constructs a y 2 function 
of the form 


X 
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where (pj eco ); are the scalar transverse momenta of the reconstructed objects, such as jets and 
photons that pass the above mentioned identification criteria, the (cy T ), are the expected reso¬ 
lutions in each object, the cr Em ; ss are the resolution of the Ey ISS projection along the x-axis and 

x,y 

the y-axis, and the (pr); are the free parameters allowed to vary in the minimization of the y 2 
function. The E™ ss terms are functions of the free parameters p x ,y, 
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In events with no genuine Ey ISS , the mismeasured quantities are re-distributed back into the 
particle momenta, to minimize the y 2 value. Events are rejected if the minimized Ey ISS (Ey ISS ) 
is less than 45 GeV and the chi-square probability is larger than 10 2 . 

To further suppress multijet backgrounds, events are not considered if the scalar sum of the 
transverse momenta of the identified jets in the event ( Hj ) is greater than 100 GeV. An addi¬ 
tional requirement is applied on the angle (a) between the beam direction and the major axis of 
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the supercluster l28l in order to reject non-prompt photons that have showers elongated along 
the beam line. 

Finally, the transverse mass, m'j 1 = V2£y£y IS5 [l — cos Ac/hy, Ey ISS )], formed by the photon 
candidate, E™ 1SS , and their opening angle, is required to be greater than 100 GeV. Photons from 
the continuum Z 7 background have a harder spectrum than the photons resulting from the 
Higgs decay in the SUSY benchmark models considered. To further reduce the continuum 
Z 7 background and for models with »?h = 125 GeV a cut of Ej < 60 GeV is applied. For 
higher masses the cut is optimized depending on each mass hypothesis going from 60 GeV up 
to 200 GeV for Mh = 400 GeV. 

The list of selection criteria used in the model independent and the SUSY benchmark model 
analyses are given in Table [TJ together with the cumulative efficiencies relative to the preselec¬ 
tion for signal and background processes. 

Table 1: Summary of ggH selection for both the quasi model-independent analysis and the anal¬ 
ysis with the SUSY benchmark model with the cumulative efficiencies of the selection require¬ 
ments relative to the preselection for Z 7 —> uvj, 7 +jet and for a signal in a SUSY benchmark 
model with ggH production of a Higgs boson with mass 125 GeV decaying into a neutralino of 


mass 120 GeV and a photon. 

Selection requirements 

Model-independent 

SUSY benchmark model 

Z 7 —> 1 / 1/7 

7 +jet 

Z 7 —> 1 / 1/7 

7 +jet 

m ~0 = 120 GeV 

A 1 

Number of jets <2 

0.909 

0.769 

— 

— 

— 

A^)( 7 ,jet) < 2.5 radians 

0.834 

0.262 

— 

— 

— 

Transverse mass > 100 GeV 

— 

— 

0.867 

0.292 

0.829 

H T < 100 GeV 

— 

— 

0.785 

0.188 

0.804 

Ej ]SS > 45 GeV 

— 

— 

0.761 

0.071 

0.743 

Prob(^ 2 ) < 10 ~ 3 

— 

— 

0.626 

0.033 

0.467 

E™ 1SS significance > 20 

— 

— 

0.440 

0.001 

0.195 

a > 1.2 

— 

— 

0.390 

0.001 

0.165 

Ej < 60 GeV 

— 

— 

0.074 

0.0002 

0.106 


4.2 Event selection in the ZH channel 

The leptonic decays of the Z boson, consisting of two oppositely charged same-flavor high-px 
isolated leptons (e ! e , p~p ), are used to tag the Higgs boson candidate events. Large missing 
transverse energy from the undetectable particles, at least one isolated high-Ej photon, and 
little or moderate jet activity are required to select the signal events. 

The details of the lepton candidate selection and missing transverse energy reconstruction are 
given in Ref. |32| . In addition, photon requirements based on a multivariate selection discussed 
in Refs. H28l[33 1 have been used. The kinematic selection requires two leptons with pj > 20 GeV 
and one photon with Ej > 20 GeV. Furthermore, the dilepton mass must be compatible with 
that of a Z boson within 15 GeV of the pole mass. 

To reduce the background from WZ events, events are removed if an additional loosely identi¬ 
fied lepton is reconstructed with px > 10 GeV. To reject most of the top-quark background, an 
event is rejected if it passes the b-tagging selection (anti b-tagging) or if there is a selected jet 
with px larger than 30 GeV (jet veto). The b-tagging selection is based on the presence of a muon 
in the event from the semileptonic decay of a bottom-quark, and on the impact parameters of 
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the constituent tracks in jets containing decays of bottom-quarks Il34t . The set of b-tagging veto 
criteria retain about 95% of the light-quark jets, while rejecting about 70% of the b-jets. 

The signal topology is characterized by a Z (££) system with large transverse momentum bal¬ 
anced in the transverse plane by a E™ 1SS + Ej system from the Higgs boson decay To re¬ 
ject background from Z 7 and Z+jets events with misreconstructed E™ ISS the azimuthal an- 

gmiss_|_g? 

gle A<^ gmissigT is required to be greater than 2.7 radians, the variable p T 1 T — pj | / pj is 
required to be smaller than 0.5, and the azimuthal angle between the two leptons A tpu is re¬ 
quired to be smaller than 2.25 radians. Finally, pj is required to be larger than 60 GeV, and 
£™ 1SS is required to be larger than 60 GeV. A summary of the selection for the analysis is shown 
in Table 12 

The signal-to-background fraction depends on the \ >f r \, the pseudorapidity of the photon, with 
greater discrimination at lower values. To exploit this effect and improve sensitivity, the se¬ 
lected events are subdivided according to whether the photon is reconstructed in the barrel or 
endcap regions, as explained in Section [Z2j 


Table 2: Summary of ZH selection. 


Variable 

Selection 

Leptons 

2 leptons, px > 20 GeV 

Photons 

1 photon, Ej > 20 GeV 

| m u - m z \ 

<15 GeV 

Anti b-tagging 

applied 

Jet veto 

0 jets with Pj > 30 GeV 

5+ s +Et 

>2.7 radians 


<0.50 

A (pu 

<2.25 radians 

Pt 

>60 GeV 

£miss 

>60 GeV 


5 Background estimation 

The background estimation techniques and the composition of all backgrounds in the search 
with the ggH and ZH signatures are discussed below. The yield for the irreducible background 
from H —» Z 7 —> i'i '7 is negligible and is therefore ignored in the analysis. 

5.1 Background estimation in the ggH channel 

The dominant background for the 7 +E™ ISS signal in the ggH channel is the process 7 +jet. Other 
SM backgrounds include Z 7 —» 1 / 1 / 7 , W 7 , W —> et/, W —» pv, W —> tv, multijet, and dipho¬ 
ton events. Background events that do not arise from pp collisions are also considered in the 
analysis. These backgrounds can be categorized broadly into three categories, as described 
below. 

5.1.1 Background estimates from simulation 

The 7 +jet process surviving the various E™ ISS selection requirements is one of the most sig¬ 
nificant backgrounds in this analysis due to the presence of an isolated photon and its large 
production cross section. The MC normalization of this background is corrected using control 
samples in data for two event classes, events without jets and those with one or more jets. The 
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control samples in data are obtained using events collected with a prescaled single-photon trig¬ 
ger and with the E™ ISS requirement reversed to ensure orthogonality to the signal phase space, 
where the contamination from other processes is minimal. Multiplicative correction factors (C) 
are obtained after normalizing the event yield in the simulation to match the data in the control 
region, separately for events with no jets (C = 1.7) and one or more jets (C = 1.1). These cor¬ 
rection factors are used to normalize the simulated 7 +jet event yield in the signal region. An 
uncertainty of 16% is obtained for these correction factors based on the difference between the 
corrected and uncorrected simulation and the relative fraction of no jet events (about 10 % of the 
events in the control region) and one or more jet events. The background processes Z 7 —t ££7 
and W —> }iv contribute only a small fraction of the total background prediction, due to the 
lepton veto applied at the preselection stage, and are modeled using simulated samples. 

To take into account differences between data and simulation due to imperfect MC modeling, 
various scale factors (SF) are applied to correct the estimates from simulation. These SFs are 
defined by the ratio of the efficiency in data to the efficiency in simulation for a given selection. 
The SF for photon reconstruction and identification is estimated from Z —> e ! e~ decays 8351 
and is consistent with unity 

5.1.2 Background estimates from data 

The contamination from jets misidentified as photons (jet —t 7 ) is estimated in a data con¬ 
trol sample enriched with multijet events defined by Ej ,ss < 40 GeV. This sample is used to 
measure the ratio of the number of candidates that pass the photon identification criteria to 
those failing the isolation requirements. The numerator of this ratio is further corrected for the 
photon contamination due to direct photon production using an isolation side band. The cor¬ 
rected ratio is applied to data events which pass the denominator selection and all other event 
requirements in the signal region. 

The systematic uncertainty of this method is dominated by the choice of the isolation sideband, 
and is estimated to be 35% by changing the isolation criteria in the sideband region definition. 
The other sources of systematic uncertainty are determined by changing the E™ ISS selection for 
the control region, and the loose identification requirements on the photons, all of which are 
found to be comparatively small. 

Events with single electrons misidentified as photons (electron —t 7 ) are another major source 
of background. This background is estimated with a tag-and-probe method using Z —> e 1 e 
events l36l . The efficiency to identify electrons (e 7e ) is estimated in the Z boson peak mass 
window of 60-120 GeV. The inefficiency (1 — e 7e ) is found to be 2.31 ± 0.03%. The ratio (1 — 
e 7e )/ c 7c , which represents the electron misidentification rate, is applied to a sample where can¬ 
didates are required to have hits in the pixel detector, and is used to estimate the contamination 
in the signal region. The misidentification rate is found to be dependent on the number of ver¬ 
tices reconstructed in the event and the number of tracks associated to the selected primary 
vertex. The difference in the final yields taking this dependence into account or neglecting it, 
using the inclusive measurement of e 7e , is less than 5%. 

5.1.3 Non-collision background estimates from data 

The search is susceptible to contamination from non-collision backgrounds, which arise from 
cosmic ray interactions, spurious signals in the ECAL, and accelerator-induced secondary par¬ 
ticles. The distribution of arrival-times of photons from these backgrounds is different to that of 
prompt photons produced in hard scattering. To quantify the contamination from these back¬ 
grounds a fit is performed to the candidate-time distribution using background distributions 
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6 Summary of systematic uncertainties 


from the data Il29l . The contamination due to out-of-time background contributions is found to 
be less than one percent of the total background and is therefore not included in the final event 
yield. 

5.2 Background estimation in ZH channel 

Processes that contribute significantly to the SM expectation in the ZH channel are listed below. 


5.2.1 Non-resonant dilepton backgrounds 

The contributions from W + W~, top-quark, W + jets, and Z/ 7 * —> t + t~ processes are esti¬ 
mated by exploiting the lepton flavor symmetry in the final states of these processes |37| . The 
branching fraction to the e ± fi T final state is twice that of the e ' e or u 1 }i final states. There¬ 
fore, the e^f-i control region is used to extrapolate these backgrounds to the e ' e and j-i 1 }i 
channels. The method considers differences between the electron and muon identification effi¬ 
ciencies. The data driven estimates agree well with the number of background events expected 
when applying the same method to simulation. The small difference between the prediction 
and the obtained value using simulated events is taken as a systematic uncertainty. 

The limited number of simulated events is also considered as part of the systematic uncertainty. 
In summary, the total systematic uncertainty is about 75%. Only two events were selected in 
the data control region. 


5.2.2 Resonant background with three leptons in the final state 

The WZ —> tv It process dominates the resonant backgrounds with three leptons in the final 
state. The electron —> 7 misidentification rate is measured in Z —)• e 1 e events by compar¬ 
ing the ratios of electron-electron versus electron-photon pairs in data and in simulation, as 
described in Section 5.1.2 The average misidentification rate is 1-2% with the larger values at 
higher |^| 7 . 


5.2.3 Resonant background with two leptons in the final state 

The WZ —> tvtt process with failure to identify the lepton from W boson decays and the ZZ —> 
2t2v process dominate these types of events. The jet —» 7 misidentification rate is measured 
in a sample containing a muon and a photon. This sample is expected to be dominated by jets 
misidentified as photons, with some contamination from W/ Z 7 events, which are subtracted 
in the study using the simulated prediction. The misidentification rate is similar to the obtained 
values in the ggH channel. 


5.2.4 Resonant background with no genuine missing transverse energy 

The background from Z 7 or Z + jets events is predicted by the simulation to be about 15% 
of the total background. Several data regions are studied to verify that the background is 
estimated correctly. A good agreement between data and simulation is found in all cases. A 
50% uncertainty, the statistical uncertainty of the control region, is taken for these backgrounds 
estimated from simulation. 


6 Summary of systematic uncertainties 

Systematic uncertainties in the background estimates from control samples in data are de¬ 
scribed in Section [5l A summary of the systematic uncertainties considered in each channel 
are listed in Tables|3]and[U 
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A common source of systematic uncertainty is associated with the measurement of the inte¬ 
grated luminosity, determined to 2.6% lf38l . The uncertainties in the normalization of signal 
and simulation-based backgrounds are obtained by varying the renormalization and factoriza¬ 
tion scales, and the parton distribution functions H26ll39ti43| . 

Because the model-independent and model-specific selections differ significantly in the ggH 
channel, the systematic uncertainties are evaluated separately for each selection. The photon 
energy scale uncertainty |28l of about 1% affects the signal and background predictions by 4% 
for the model specific selection and by 0.5% for the model-independent selection. Similarly, the 
jet energy scale uncertainty affects the signal and background predictions by 2-5% depending 
on the process and selection. After changing the photon or jet energy scales, the £j ,ss is also 
recomputed. In addition, the systematic uncertainty associated with the jet energy resolution 
and unclustered energy scale are propagated to the Ey ISS computation, and affect the signal 
and background predictions by 2-4%. As described in the previous section, a 16% uncertainty 
is applied to the 7 +jet normalization due to the difference in the jet multiplicity distribution 
between the data and background predictions in the 7 +jet control region. The uncertainty due 
to the pileup modeling is found to be 1 %, and is estimated by shifting the central value of the 
total inelastic cross section within its uncertainty. 


Table 3: Summary of all relative systematic uncertainties in percent for the signal and back¬ 
ground estimates for the Higgs model (model-independent in parenthesis) selection in the ggH 
analysis. 


Source 

Signal 

Jet—* 7 

Electron—> 7 

7 + jet 

Z++7 

W 7 

PDF 

10 ( 0 ) 

— 

— 

— 

4(4) 

4(4) 

Integrated luminosity 

2.6 ( 2 . 6 ) 

— 

— 

2.6 ( 2 . 6 ) 

2.6 ( 2 . 6 ) 

2.6 ( 2 . 6 ) 

Photon efficiency 

3(3) 

— 

— 

3(3) 

3(3) 

3(3) 

Photon energy scale ± 1 % 

4 (0.5) 

— 

— 

4 (0.5) 

4 (0.5) 

4 (0.5) 

Ey ISS energy scale 

4(2) 

— 

— 

4(2) 

4(2) 

4(2) 

Jet energy scale 

3(2) 

— 

— 

5(5) 

3(2) 

3(2) 

Pileup 

1 ( 1 ) 

— 

— 

1 ( 1 ) 

1 ( 1 ) 

1 ( 1 ) 

Z 1 / 1/7 normalization 

— 

— 

— 

— 

3(3) 

— 

7 + jet normalization 

— 

— 

— 

16 (16) 

— 

— 

W 7 normalization 

— 

— 

— 

— 

— 

3(3) 

Jet —> 7 

— 

35 (35) 

— 

— 

— 

— 

Electron —> 7 

— 

— 

6 ( 6 ) 

— 

— 

— 


In the ZH channel, lepton-reconstruction and identification scale factors are measured using a 
control sample of Z/ 7 * —> i + i~ events in the Z peak region |[36l . The associated uncertainty is 
about 2% per lepton. The photon identification uncertainty is taken to be 3% l33l . The effect of 
uncertainties in jet-energy scale and E™ ISS on the analysis is also considered. The uncertainty 
in the b-tagging efficiency is estimated to be about 0.7% comparing inclusive Z/ 7 * —» £ + £~ 
samples in data and simulation. The total uncertainty in the background estimates in the signal 
region is 36%, which is dominated by the statistical uncertainties in the data control samples 
from which they are derived. 

The impact of the systematic uncertainties in the ggH channel is relatively important: the sen¬ 
sitivity is increased by about 50% if all the systematic uncertainties are removed, where the 
normalization uncertainties on the 7 + jet and jet —» 7 background processes dominate. The 
ZH channel is limited by the statistical uncertainty, and the effect of the systematic uncertainties 
reduce the sensitivity by less than 10 %. 
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7 Results 


Table 4: Summary of relative systematic uncertainties in percent for the signal and background 
estimates in the ZH analysis. 


Source 

ZH 

Z 7 or Z + jets 

WZ 

zz 

WW + top-quark 

Integrated luminosity 

2.6 

— 

2.6 

2.6 

— 

Lepton efficiency 

3.6 

— 

3.6 

3.6 

— 

Photon efficiency 

3.0 

— 

— 

— 

— 

Momentum resolution 

0.5 

— 

1 

1 

— 

Ey ISS energy scale 

0.5 

— 

0.6 

0.1 

— 

Jet energy scale 

2 

— 

4 

4 

— 

b-tagging 

0.7 

— 

0.7 

0.7 

— 

Underlying event 

3 

— 

— 

— 

— 

PDF 

7.1 

— 

6.3 

7.7 

— 

Renorm, and factor, scales 

7.0 

— 

10.7 

6.5 

— 

Z/ 7 * —> l + £~ normalization 

— 

50 

— 

— 

— 

Non-resonant dilepton bkg. norm. 

— 

— 

— 

— 

70 

Jet —> 7 

— 

— 

30 

30 

— 

Electron —> 7 

— 

— 

10 

10 

— 

Amount of simulated events 

3.5 

60 

10 

30 

40 


Correlations between systematic uncertainties in the two channels are taken into account. In 
particular, the main sources of correlated systematic uncertainties are those in the experimental 
measurements such as the integrated luminosity, photon identification, the jet energy scale, 
and missing transverse energy resolution. All other systematic uncertainties are uncorrelated 
between them given the different signal and background processes. 


7 Results 

The results from the two searches and their combination are reported in this section. In the 
absence of deviations from the standard model predictions, the modified frequentist method, 
CL S Ii44l - l46l . is used to define the exclusion limits. 

7.1 Model-independent results in the ggH channel 

Because of the variety of possible BSM signals that could contribute to this final state, the results 
are presented for a signal with the model-independent selection described in Section |4] The 
total number of observed events and the estimated SM backgrounds are summarized in Table[5j 

rn y£miss 

and found to be compatible within their uncertainties. Figure 2 shows the nij T and Ey ISS 
distributions for the model-independent selection. 

Figure [3] shows the observed and expected model-independent 95% CL upper limits for the 

y£ miss 

ggH analysis on the product of cross section, acceptance, and efficiency for mj T > 100 GeV, 
as a function of E™ ISS threshold. 

7.2 Model-specific results in the ggH channel 

Imposing the model-specific selection described in Section [4] for the ggH channel, 1296 events 
are selected in data with a total estimated background of 1232 ± 188. The yields for this se¬ 
lection are shown in Table [o] and estimated for Higgs boson decays (H —> Gy f |, x'f —> G 7 ) 
assuming the ggH production rate for SM Higgs bosons and a 100% branching fraction for this 
decay. Figure [4] shows the transverse energy distribution of photons for data, the background 
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Table 5: Observed yields and background estimates at 8 TeV in the ggH channel after the model- 
independent selection. Statistical and systematic uncertainties are shown. 


Process 

Event yields 

7 + jets 

(313 ± 50) x 10 3 

jet —> 7 

(910 ± 320) x 10 2 

e —>■ 7 

10350 db 620 

W (—> £v) + 7 

2239 ± 111 

Z(-> vv) + 7 

2050 ± 102 

Other 

1809 ± 91 

Total background 

(420 ± 82) x 10 3 

Data 

442 x 10 J 




CMS 


7.4 fb 1 (8 TeV) 
y + Jets 
Jet -> y 
Electron -> y 
Z y 

W y ->l v y 
W —» |i(x) v & y Y 

Data 

ggH (125,120) 


CMS 


7.4 fb' 1 (8 TeV) 


Y + Jets 
Jet -> y 
E lectron -» y 
Z y 

W y —> I v y 
W —» fi('u) v & y Y 

Data 

ggH (125,120) 


icr 1 

d. 2 

m 1 

S’ 0 

“_1 

ra 

ra¬ 


ce-2 


100 


150 


200 


250 300 

i-miss 


[GeV] 


m 10 

°io 7 

O 

^io 6 

|io 5 

£l 0 * 

LU 

10 3 

10 2 

10 


50 100 150 200 250 300 

m T [GeV] 


^£miss 

Figure 2: The nij T and £y ISS distributions for data, background estimates, and signal after 
the model-independent selection for the ggH channel. The bottom panels in each plot show the 
ratio of (data - background)/background and the gray band includes both the statistical and 
systematic uncertainties in the background prediction. The signal is shown for mu = 125 GeV 
and m ~o = 120 GeV and a 100% branching fraction. 
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7 Results 


CMS 7.4 fb' 1 (8 TeV) 

^ L I I I I I I I I I I I I I I I I I I I I I I I I I I I - 

^ - M T > 100 GeV 

<a - 10 3 - - 95% CL observed ~ 

D : 


10 2 


10 


Figure 3: The expected and observed 95% CL upper limit on the product of cross section, ac- 

-,£miss 

ceptance, and efficiency (cr(pp —» 7 + E™ lss ) Ae) for nij T > 100 GeV, as function of the Ey ISS 
threshold for the ggH channel. 

estimates, and signal after the model-dependent selection, except the upper selection on the 
photon, for the ggH channel. 

7.3 Results in the ZH channel and combinations 

A total of four events are selected with the search in ZH. The background yield is estimated 
to 4.1 ± 1.8. The numbers of observed and expected events are shown in Table [7j The signal 
model assumes a SM ZH production rate and a 100% branching fraction to undetectable parti¬ 
cles and one or two photons. The expected signal yield is larger for cases where m ~o is smaller 

than »zh/ 2 since there are two photons in the final state (H —> j y ( l , y ( l ) —> 77 GG), and as a result 
the sensitivity improves for smaller masses. Good agreement between the data and the back¬ 
ground prediction is observed. The transverse mass, mj Ej = \flpjPj T + ^ T [1 — cos(A<^ gmissigr)], 
and \t ] 7 1 distributions discriminate signal and background and are shown in Figure [5]at the fi¬ 
nal step of the selection. 

The 95% CL upper limits are extracted from counting experiments in three categories: the 
model-specific selection in the ggH channel, and photons identified in the barrel and the end- 
cap calorimeters for the ZH channel. Results are combined using a binned-likelihood method. 

The 95% CL upper limits on (a B) / <Jsm, where is the cross section for the SM Higgs boson, 

are evaluated for different mass values of x® ranging from 1 GeV to 120 GeV for the individual 
searches and their combination and are shown in Figure | 6 j The upper limits for m-o < mu/2 
are not shown for the ggH channel because the sensitivity is very low due to the combination 
kinematic properties and the corresponding selection; in particular the E™ ISS and photon pj 
values tend to be outside the selected ranges. A 95% CL upper limit on the branching fraction 
of 10% is set for a neutralino mass of 95 GeV. 

Expected and observed limits are also shown for the decay of possible heavier scalar Higgs 
bosons as a function of the Higgs boson mass in Figure Wj The requirement on Ey used in 
the ggH channel is removed. A lower threshold on Ey is added, optimized to maximize the 


My > 100 GeV 

-95% CL observed ~ 

. 95% CL expected ; 

^ + 1o expected 
| ± 2 g expected 



6s ?S 6s 9s ro S ’7s 7 SS ’3S ’<S 1 SS 7 6S 7 ?S 7 6S 7 9S 

Lower Ey lss Threshold [GeV] 
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Table 6 : Observed yields, background estimates, and signal predictions at 8 TeV in the ggH 
channel for different values of the m^o and for different ct~o of the y) 1 . These correspond to 
£>(H —>■ undetectable + 7 ) =100%, assuming the SM cross section at the given ni\\ hypothesis. 
The combination of statistical and systematic uncertainties is shown for the yields. 


Process 


Event yields 


ggH(mn = 125 GeV, rn~o = 65GeV) 
ggH(mn = 125 GeV, rn-o = 95 GeV) 
ggH(m H = 125 GeV, = 120 GeV) 


653 ± 77 
1158 ± 137 
2935 ± 349 


ggH (m H = 125 GeV, m^ = 95 GeV) ct~o = 100 mm 983 dh 116 

ggH(mn = 125 GeV, m^o = 95 GeV) ct~o = 1000 mm 463 ± 55 

ggH(mn = 125 GeV, m^o = 95 GeV) ct~o = 10000 mm 83 ± 10 

ggH(m H = 150 GeV, rn^ = 120 GeV) ' 4160 ± 491 

ggH(m H = 200 GeV, m^> = 170 GeV) 5963 ± 704 

ggH(m H = 300 GeV, = 270 GeV) 5152 ± 608 

ggH(??z H = 400 GeV, = 370 GeV) 4057 ± 479 


7 + jets 
jet -4 7 
e —> 7 

W(—* tv) + 7 
Z(—> vv) + 7 

Other 


179 ± 28 
269 ± 94 
355 ± 28 
154 ± 15 
182 ± 13 
91 ± 10 


Total background 


1232 ± 188 

Data 



1296 


CMS 

7.4 fb' 1 (8 TeV) 

10 6 


1 1 ■ y + Jets 

Jet -> y 



Electron -4 y 
Z 7 

|W y ->l v y 
|W 4 |x(x) v & y y 

• Data 

ggH (125,120) 



1 



~ f- . ^ . . j 


1 i ■ t it I it i it it 1 it i >r it 1 it i r -i I i 


a 50 100 150 200 250 300 350 400 450 500 

[GeV] 


Figure 4: The transverse energy distribution of photons for data, the background estimates, 
and signal after the model-dependent selection (except the upper selection on the photon) for 
the ggH channel. The bottom panel shows the ratio of (data - background) /background and 
the gray band includes both the statistical and systematic uncertainties in the background pre¬ 
diction. The signal is shown for niy\ = 125 GeV and rn-o = 120 GeV. 
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7 Results 


Table 7: Observed yields, background estimates, and signal predictions at 8 TeV in the ZH 
channel for different values of the and for different ct~o of the x^. The signal predictions 
correspond to 15 (H —> undetectable + 7 ) =100% assuming the SM ZH cross section at the 
given hypothesis. The combination of statistical and systematic uncertainties is shown for 
the yields. 


Process 

ZH(7?zh = 125 GeV, 777-0 
ZH(/?z h = 125 GeV, 7/7-0 
ZH(m H = 125 GeV, /77-0 
ZH(m H = 125 GeV, /77-0 
ZH(m H = 125 GeV, /77-0 
ZH(m H = 125 GeV, /77-0 
ZH(/ 7 z h = 125 GeV, 777-0 
ZH(77Zh = 125 GeV, 777-0 
ZH(77Zh = 125 GeV, 777-0 
ZH( 77 z h = 125 GeV, 777-0 
ZH( 777 h = 200 GeV, 777^0 
ZH( 77 z h = 300 GeV, ?77^o 
ZH( 77 z h = 400 GeV, 777^0 
Z 7 + Z + jets 

wz 

zz 


= 1 GeV) 

= 10 GeV) 
= 30 GeV) 
= 60 GeV) 
= 65 GeV) 
= 95 GeV) 
= 120 GeV) 


= 95 GeV) c i~o = 100 mm 
= 95 GeV) ct~o = 1000 mm 
= 95 GeV) ci~o = 10000 mm 
= 170 GeV) 

= 270 GeV) 

= 370 GeV) 


WW + top-quark 
Total background 
Data 


Event yields 

69.2 ± 8.4 

68.6 ± 8.4 

53.5 ± 6.5 

47.7 ± 5.8 
40.0 ± 4.9 

40.3 ± 4.9 
39.0 ± 4.8 

39.3 ± 4.8 

17.6 ± 2.2 

2.6 ± 0.3 

13.1 ± 1.6 
3.5 ± 0.4 
1.2 ± 0.1 

0.6 ± 0.4 
1.2 ± 0.3 
0.3 ± 0.1 
2.0 ± 1.7 

4.1 ± 1.8 
4 
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Figure 5: Distributions in signal where Wh = 125 GeV and m-o = 95 GeV, backgrounds and 

£££miss 

data for m T T (left) and | ;/ 7 1 (right) after applying all requirements. The uncertainty band 
for the backgrounds includes both statistical and systematic uncertainties. The signal model 
assumes a SM ZH production rate for a Higgs boson with m h = 125 GeV and a 10% branching 
fraction. 


sensitivity for each mass hyposthesis. A combination of the two channels is not performed 
because the assumption of a common SM Higgs boson cross section is not justified. 

As discussed in the introduction, some BSM models predict Xi neutralinos with sizable life¬ 
times. The performance of the searches has been evaluated for finite lifetimes without modify¬ 
ing the analysis strategy The expected and observed limits are shown in Figure [8] as function 
of ct~o. The results are shown for /»h = 125 GeV and m-o = 95 GeV. As seen in Tables |6j and pj 
the selection efficiency is roughly constant for values of ci~o less than 10 cm, and drops rapidly 
for larger values. The default timing criteria applied in the ECAL energy reconstruction are 
the cause for the decrease in the efficiency In particular, there is a requirement of a maximum 
of 3 ns on the photon arrival time relative to the nominal time-of-flight for prompt photons. 
The delayed arrival time of the photon can be caused by a kink in the trajectory or by a lower 
velocity of the neutralino. 


8 Summary 

A search is presented for exotic decays of a Higgs boson into undetectable particles and one or 
two isolated photons in pp collisions at a center-of-mass energy of 8 TeV. The data correspond 
to an integrated luminosity of up to 19.4 fb 1 collected with the CMS detector at the LHC. 
Higgs bosons produced in gluon-gluon fusion or in association with a Z boson are investigated. 
Models including Higgs boson decays into a gravitino and a neutralino or a pair of neutralinos, 
followed by the neutralino decay to a gravitino and a photon, are tested. The measurements for 
the selected events in data are consistent with the background only hypothesis, and the results 
are interpreted as limits on the product of cross sections and branching fractions. Assuming a 
standard model Higgs production cross-section, a 95% CL upper limit is set on the branching 
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8 Summary 



Figure 6 : Expected and observed 95% CL upper limits on ctB/ctsm for in H = 125 GeV as a 
function of m^o assuming the SM Higgs boson cross sections, for the ZH and ggH channels and 

their combination, with B = B (H —» XiX^_) -A G + y ) 2 for m^o < m.n/2 and B = B (H —> 
^?G) B -A G + 7 ) for m^o > m n /2. 

fraction of a 125 GeV Higgs boson decaying into undetectable particles and one or two isolated 
photons as a function of the neutralino mass. For neutralino masses from 1 to 120 GeV an 
upper limit in the range of 7 to 13% is obtained. Further results are given as a function of the 
neutralino lifetime, and also for a range of Higgs boson masses. 
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Figure 7: Expected and observed 95% CL upper limits on h B as a function of the Higgs 
boson mass with = m\\ — 30 GeV in ggH channel (left) and in the ZH channel (right). 



Figure 8: Expected and observed 95% CL upper limits on Uh B as a function of ct~o for zh h = 
125GeV and m~ o = 95 GeV, where B = £?(H —> x?G) B(x ® —> G + 7 ). 
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